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Abstract 

We start from a Hamiltonian describing non-interacting electrons and add bosons 
to the model, with a Jaynes-Cummings-like interaction between the bosons and 
electrons. Because of the specific form of the interaction the model can be solved 
exactly. In the ground state, part of the electrons form bound pairs with oppo- 
site momentum and spin. The model also shows a gap in the kinetic energy of 
the electrons, but not in the spectrum of the full Hamiltonian. This gap is not 
of a mean-field nature, but is due to the Pauli exclusion principle. 



1. Introduction 

In the theory of strongly correlated particles a model known for a long time 
is the mean-field model with Hamiltonian 

k^l k 

This is the mean-field BCS model [H, 0, S @1 • Despite the simplicity of this 
model many interesting results about strongly coupled particles can be derived 
from it. Recently, however, there has been a lot of interest in other models 
for the creation of paired electrons, motivated in large part by the problem of 
finding a theory for high T c superconductivity 

Among the most useful models for stron gly correlated systems are the Hub- 



bard model and its generalizations [8|, lil. [ToTTlll. These models have been used 
to describe not only superconductivity [1 13L 14 1 but also ferro mag netism [IF 



metal-insulator transitions [16J and other properties of solids [17|, uM an d of 
atoms in optical lattices (l9j . 

In this paper we start from a Hamiltonian describing the motion of non- 
interacting electrons on a lattice. We then add to the model bosons, which 



'Corresponding author 
Email addresses: tobias.verhulstSua.ac.be (T. Verhulst), jan.naudtsOua.ac.be (J. 
Naudts) 

1 Research Assistant of the Research Foundation - Flanders (f WO - Vlaanderen) 



Preprint submitted to Elsevier 



May 10, 2010 



interact with the electrons in a way similar to the boson-fermion interaction in 
the Jaynes-Cummings model [20| . The complete Hamiltonian is then similar to 



the one used to describe bipolaron interactions |2l|, |22j, |23[. However, in our 
model this interaction term does not describe the interaction between pairs of 
electrons, but rather provides a process which creates and destroys pairs. A 
short account of the present findings has already been given in 24 1 . 

In models like £T|) the binding energy of all pairs is the same. It has recently 
been suggested [25], |26| that, due to the Pauli exclusion principle, this energy 
should not be a constant. In our model we find that, indeed, the binding energy 
depends on the momenta k, — k of the paired electrons. 

Overview of the paper. In section 2 we give a detailed description of the 
different terms in the Hamiltonian, both their mathematical properties and 
physical interpretations. In section 3 the exact solution of the model for a 
fixed value of k, — k is presented. In section 4 the many particle ground state 
is constructed. In section 5 the elementary excitations and the gap in the 
energy are discussed, followed by the main conclusions in section 6. A list of all 
eigenstates and eigenvalues of the model can be found in the appendix. 



2. The model 

The Hamiltonian we consider consist of three parts: one describing the elec- 
trons, one describing the bosons and one describing an interaction between 
them. The fermion part is given by 

Hp = -hJ2 Akfll, A„ + Bt >(T B k<a Bl_ k<T B_ K „ + Wk-fr (2) 

k,cr k.cr k 

The first term comes from the kinetic energy T — — ■ tij b\ a bj^ which is 
written in diagonal form using Bogoliubov's quasi-particles. The quasi-particle 
operators are defined by B\^ a — J^. Uk,i&i,<7 with v^i the matrix coefficients of 
the diagonalizing unitary transformation. These operators satisfy the canonical 
anti-commutation relations. The Ak are the eigenvalues of the kinetic energy. 
On a cubic N x N x N lattice they are given by 

3 

A k = 2 ^2 cos(27rfc a ) (3) 

a=l 

where k a = 2 "~^ +1 for n from to N - 1. 

The second part of Hp describes an interaction between electrons with equal 

spin. If the parameter 77 is positive, this term raises the energy of states where 

electrons with opposite momenta have equal spins. 

The operator J£ in the last term of fl5J is defined as J£ — > J^] with 
= B_ k ^B^ ^B-k t ±Bkrf and j£ = (J±)^. and together generate a 

representation of su(2) and can be interpreted as some kind of spin operators. 
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An expression for J k in counting operators nk,o- = a B] ia can be found: 

•K = ^"-k,t(! - n_k,t)«-k,4.(l - n-k,;) - ^"-k,^ 1 - "-k,;) n -k,t( 1 ~ n k,t)- ( 4 ) 

It is then easy to see that J k has eigenvalues —4,0, and 5, as one would expect. 
This part of the Hamiltonian describes the interaction with a constant external 
field. 

Notice that one could add some further terms quadratic in the -Beer's. For 
instance, we do not consider a spin flip part in the Hamiltonian, or terms chang- 
ing the total momentum. We could add a term proportional to ^ k a B-^ a . 
This would make the calculations more complicated, although the model remains 
integrable. However, such a term does not change anything fundamentally. So 
we neglect it here. 

The bosonic part of the Hamiltonian describes a collection of harmonic oscil- 
lators, one for each fermionic degree of freedom. The creation and annihilation 
operators are denoted Ck, r and c t T , with r =| or =4,. For simplicity of notation 
we use the same labels as for the fermionic operators. These oscillators could 
for instance be lattice vibrations with two possible polarizations. The bosonic 
Hamiltonian reads 

h b = ^ ^ k ( c k,T ck <t + c k,4 c k,i) = X! ^ k X! mk ' T ' ( 5 ) 

k ' k r=t,l 

with /ik = |U_k > and with mk. T = c t r Ck. T - 

By introducing new operators = Ck,f c_ k , this bosonic Hamiltonian can 
be rewritten as 

H B = ^ ^ k ( a k°k + Oka k ) _ ^A* kTOk ,t m -k4- (6) 

k k 

While the operators Ck,r and c k satisfy the canonical commutation relations 

for bosons, the ak and a k operators do not. It still holds that [oik, ai] = 0, but 
the commutation between the new creation and annihilation operators is more 
complicated: 

[a k , a[] = <5 k ,i (-Ck,t4:,t + c -*>l c -k,l) = ( TOk ,t ~ TO -k,i) • (7) 

The operators ak and a k are introduced because they are useful to express the 
interaction term (J8J), which will now be introduced. 

The interaction between the electrons and bosons is described by the Hamil- 
tonian 

^/ = ^E^( a k j k +« k j k + ) (8) 

k 

with ak, a k and as defined above. This interaction is of fourth order in the 
fermionic operators i?k,cr and , and of second order in the bosonic operators 
c k , CT and c f k cr . 
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The complete Hamiltonian H = Hf+Hb+Hi describes electrons moving on 
a lattice and interacting with each other both directly and mediated by lattice 
vibrations. Notice that the latter interaction is of the Jaynes-Cummings type 
and is expressed in the operators Jk and ctk, not in the original creation and 
annihilation operators -Bk,c and Ck,o-- 

We can compare the present model with the mean-field model (HJ by identi- 
fying JjT with cr^" <r~ k . Here, there are only interactions between particles when 
they have the same or opposite momentum ±k, while in (TTJ) there is a mean- 
field interaction between particles with all possible momenta. In addition, the 
interaction term of our model is of second order in the Ok 



Thus, integrating out the bosons gives a term of fourth order in the <7k 
instead of one of second order. 

3. The exact solution of the model 

Because there is no interaction between particles with momenta k respec- 
tively k' unless k = ±k', we construct solutions in a subspace "Hk with only 
particles with the same or opposite momenta. Thus, the number of electrons is 
between zero and four while the number of bosons is between zero and infinity. 

We use the following basis vectors: 



where n = or 1 is the number of electrons for the given momentum and spin, 
and m counts the number of corresponding bosons. There are two important 
things to note about these basis vectors. First: most terms of the Hamiltonian 
are already diagonal in this basis. Only the fermion-boson interaction Hi is 
not. Second: on the vectors where it is not zero, the interaction Hamiltonian 
Hj acts like a Jaynes-Cummings Hamiltonian. This means that the well known 
techniques to solve the Jaynes-Cummings model can be applied. 

3.1. The diagonal parts of H 

Most terms of H are diagonal in the basis (fTOf . The eigenvalues are given 
here. For the kinetic energy of the electrons 





k 



|»k,t> m k,t! n k,4-> m k,|5 "-k,t> ™-k,t; n -k,i, TO-k,4.) 



(10) 



T = — 



(ii) 



the eigenvalues are 



7iAk(nk,t + "k,4.) - 7iA_k(ri-kt + n -^,i)- 



(12) 
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The term | J^u w k<^k nas eigenvalues 



fi(w k + u;_ k ) 

2 i n k,t(! - "-k.tJ^-kal 1 - n k,|) 

-«k,^(l - "--k,i)"-k,t( 1 - n k,t)] • (!3) 
The second term of Iff, the fermion-fermion interaction, has eigenvalues 

r/ft("k,t™-k,t + «-k,jn-k,i) . (14) 
The eigenvalues of the bosonic Hamiltonian ffg arc 

Mk (w k ,t + m-k,j.) • (15) 
All eigenstates and eigenvalues are listed in the appendix. 

3.2. The interaction Hamiltonian Hi 

The interaction Hamiltonian Hi vanishes when the number of electrons is 
not equal to two. In addition, straightforward calculation shows that it vanishes 
when the following conditions are satisfied 

1- n k,t = or ra-k,t = or m_k4 = 0, 

2- «-k,t = or nk,t = or mk,f = 0, 

3- «-k.t = or 7ik,t = or mk,j. = 0, 

4- rik.t = or n_k,t = or m_k,t = 0- 

In particular, the only basis vectors on which Hj does not vanish are the two- 
fermion vectors with vanishing momentum and vanishing spin 

|l,rak,t;0,mk,4.;0,m_k,t; l,m_k,4.) ( 16 ) 
|0, m k ,t; 1, m k:i ; 1, m_ k ,t; 0, m_ k ,+) . (17) 

The eigenstates of can now be constructed in a way similar to the solution of 
the Jaynes-Cummings model. Although it is not needed in order to obtain the 
eigenvalues and eigenvectors we assume from now on that Ak = A_k, as is for 
instance the case for a cubic lattice. This assumption simplifies the calculations 
a lot. The eigenvectors on which Hi does not vanish are denoted as follows 

|2,0, +k, +) = sin^|l,mk,t;0,mk,4.;0,m_k,t;l,TO_k,4. + 1) 

cos6»+|0,77i k ,t + 1; l,m k ,4,; l,m_k,t;0,m_ k ,4.) (18) 
|2, 0, +k, — ) = cos0 k ~|l,m k ,t + l;0,m k ,+;0,m_k,t; l,ra-k,4.) 

-sin6£|0,TOk,t;l,™k,4.; 1 > m -k,t; ' m -k,.l + 1) (19) 
|2,0,-k,+) = sin^|0,mk,t;l,"ik,4.;l,m_k,t;0,m_k,4. + l) 

cos ^|l,mk,t + l;0)"»k,4.;0,m_k,t; l,"i_k,j.) (20) 
|2, 0, — k, — ) = cos^ |0,m k ,t + 1; l,m k ,4.; I,m_ k ,t;0,m_k4) 

-sin0 k "|l,mk,t;O,mk,4.;O,m_k,t; l,ni-k4 + !)■ ( 21 ) 
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In these expressions the angles 9^ are given by 

0± = arctan ( ^ k + ^ k T + ^>-k) 2 + 4g) . (22) 

The eigenvalues are for the first and the third class, and for the second 
and the fourth class, with 

4 = e ho - 2^A k ± J yVk + ^- k ) 2 + 4£ 2 . (23) 
The energy -E^q of the harmonic oscillators equals 

e ho = TO k,t + ™k,i + ™-k,t + ™-k,j.- (24) 

Note that at least one of the harmonic oscillators is in an excited state. Thus, 
the contribution from the harmonic oscillators does not vanish. 

3. 3. The ground state for fixed k 

First we construct the ground state for a fixed value of k. Since particles do 
not interact except when they have equal or opposite momenta, the complete 
many-particle ground state can be constructed from the ground states for fixed 
k. We assume from now on that there is no external field, so u;±k = in ([2]). 
The eigenvalues for the Jaynes-Cummings-like states then become 

4 = 4o - 2^A k ± 7I&. (25) 

For simplicity let us assume now that £k > 0. If it is large enough then one of 
the is the lowest energy level. This requires that the contribution from the 
bosons is minimal, this is, equal to E^ Q = hfi^. Even then the ground state 
is still degenerate because of the symmetry between the momenta k and — k. 
Indeed, for large enough £k: the ground states are 

|2,0,k,+) = cos6£|l,l;0,0;0,0;l,0)-sin6>+|0,0;l,0;l,0;0,l) (26) 
2,0, -k, +) = cos6> k |0, l;l,0;l,0;0,0)-sin6£|l,0;0,0;0,0;l,l). (27) 

If £k is small the ground state is the unique state with four electrons and no 
bosons 

|4,0,0) = |1,0;1,0;1,0;1,0) (28) 

This eigenstate has eigenvalue £4 = — AhX^ + 2hrj. The critical value of £k at 
which the quantum phase transition occurs is 

6c = /-ik + 2(A k -77k). (29) 

The second term comes from the difference between and £4 with the same 
number of bosons while the first term stems from the fact that for the interacting 
states E\, n can not vanish. 
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3.4- The many-particle ground state 
Substituting © in gives: 



& = Mk + 2 (^2 ^ cos(27rfc Q ) - 77 J (30) 

or 

^ 008(2^) = I -f + |. (31) 

From this equation one can, at least in principle, determine at which k it be- 
comes energetically favorable to have a "pair" of electrons bound by a boson 
instead of four non-interacting electrons. 

In order to do this one needs, in addition to the assumption §5$ for Ak, 
expressions for £k and /_ik- For the interaction strength £k it is reasonable to 
assume no dependence on k. We do not specify an explicit dispersion relation 
for the bosons. 



4. Filling the band 

The number of electrons which occupy a certain (k, — k)-level varies between 
zero and four. If the interaction strength is zero, all electrons occupy levels with 
four electrons and zero bosons. However, with a nonzero interaction part of the 
electrons form pairs. Thus, for some (k, — k)-levels there will be two electrons 
and one boson instead of four electrons and no bosons. This means that, if the 
number of electrons stays the same, some of them must occupy higher k-levels. 
This effect inhibits to some extend the formation of pairs because some of the 
pairs have a higher kinetic energy. 

We describe in detail the case of a one dimensional system with a band 
filled up to the Fermi level at k = kp. We denote Xk F = Xp and ^k F = Vf- 
|4, 0,0)fc F is then the highest occupied level with energy AH(rj — Xp). If these 
electrons occupy two levels with two electrons each instead of one level with 
four electrons, then one pair has k = kp (the level |2, 0, k, — }) and the other one 
k = kp + Sk (the level |2, 0, k + 5k, —)). The pair at k = kp has the energy 

62 = hn F + 2h{ri - X F ) - %£, (32) 
and the pair at k = kp + Sk has the energy 

htJ,k F +6k + ZHv - Xk F +5k) - T>4- (33) 
Thus, the energy gained by creating two pairs from four free electrons equals 
Ae{k F ) = 2h(X kF+S k - X F + £) - h(n F + [ik F +8k) (34) 
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Around the Fermi level one can make expansions for A& and jik 



\ k = \p + ^-(\p)5k+{(8k) 2 ) (35) 

fi k =fiF + ^r(HF)5k+((6k) 2 ). (36) 
Substituting these expansions in (f34f gives 

Ae(fc F ) = h{2£, — fip) + h (^(A F ) - ^(/*f)) ^. (37) 

Now the four electrons occupying a second fc-level can be split into two pairs. 
The highest four fermion level is the one with k = kp —8k. The pairs are formed 
at k = kp — 5k and k = kp + 28k. The gain in energy is 

Ae(fc F - 8k) = h(2£ -np) + h (e^(A F ) - ^0*f )) 5k. (38) 

The same procedure where the particles from a four fermion level at k = kp—n8k 
are combined in two paired states at k = kp — n8k and k = kp + (n + l)<Jfc can 
now be repeated N times. N can be determined from the condition 

Ae(k F -N8k) = h(2Z-np) + ri ( 2(2N + 1)^-(A F ) - ^(/i F ) ) <5fc = 0. (39) 

\ dk dk J 

Some straightforward algebra then shows that the number of pairs formed is 



2^(A F ) 28k^(X F y 



2N + 1 = f x vf ^ J ' + " s . (40) 



Notice that, all other things being equal, a smaller 8k implies a larger N. Thus, 
the closer the possible values of k are to each other the more pairs are formed. 
In general, 8k is inversely proportional to the size of the system, so the fraction 
of pair-forming electrons is constant as a function of the system size. 



5. Excitations and kinetic energy gap 

The excitations with the lowest energy are the production of two fermion 
pairs at k = fc F — (N + l)8k and k = fc F + (N + 2)8k and the recombination 
of two pairs into one four fermion level at k — kp — N8k. During the creation 
of two pairs from the particles at k = kp — (N + 1)8 k, the electrons in the pair 
with the highest fc-value each get an additional momentum of (2N + 3)8k. The 
momentum of the pair with the lowest fc-value does not change. Thus, the total 
change of momentum is 

Ak tot = 2(2N + 3)8k. (41) 
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One concludes that there is a gap in the kinetic energy of the electrons. Substi- 
tuting (l40l) gives 



Notice that the second term does not depend on 5k, so even if 5k is very small 
Ak tot may still be macroscopic. An analogous reasoning holds for the reversed 
process in which two pairs are destroyed. 

6. Conclusions 

In summary, we start from a free hopping Hamiltonian for the electrons 
to develop a model with both electrons and bosons. The interaction between 
the bosons and the electrons is inspired by the interaction term in the Jaynes- 
Cummings model. 

Due to the interaction term part of the electrons become bound in pairs 
with opposite momentum and spin. In general, the many particle ground state 
contains both paired and unpaired electrons. Equation (14TJ1) shows in the one- 
dimensional case how the number of pairs formed in the ground state depends 
on the parameters of the model. 

We also demonstrate the existence of a gap in the kinetic energy spectrum. 
The gap in this model is not of a mean-field type but is a consequence of the 
Pauli exclusion principle. We also find that the binding energy of the pairs is not 
a constant but depends on the wavevector k. This is in agreement with earlier 
results [1^, [26| which suggest that, due to the Pauli exclusion principle, the 
binding energy of Cooper pairs depends on k. In the one-dimensional case, the 
total change in momentum for the elementary excitation in which two additional 
pairs are formed is given by (|4"2"j) . 

Future research is needed to study the thermodynamical properties of the 
system. Note that it is possible to write our model as a mean-field model 
similar to ((T|). This can be done by first integrating out the bosons, which 
gives a Hamiltonian of fourth order in the operators, and then introducing 
new operators and which create and annihilate two bound pairs. Using 
these new operators an effective Hamiltonian can be constructed, similar to (fT]). 
However, this new mean-field model still allows the binding energy of the pairs 
to depend on k, which is not the case in (flj. 

Appendix A. The eigenstates and eigenvalues 

Here we give an overview of all eigenstates and eigenvalues of the model. 
The notation is as follows: the first number counts the electrons and the second 
number is the total spin. The third label is the total momentum of the elec- 
trons. The quantum numbers of the bosons are omitted. In one case additional 
quantum numbers are needed. 




(42) 
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We further use the notation E\ Q for the eigenvalues of the boson Hamil- 
tonian for one k, see (l24l) . In particular, the eigenvalues of Hb have the form 

J2k E HO- 

1. The eigenvectors with vanishing electron count are of the form 

0, 0, 0} = |0, mk,t;0, m kii ;0,m_ k!t ;0,m_ k ^). (A.l) 

The energy is Sq — E^ , this is, only the harmonic oscillators contribute. 

2. With one electron there are four classes of states because the electron can 
have either momentum k or — k. They are 

l,+-,+k) = |l,mk,t;0, mk,i;0,m_k,t;0,m_k,4.), 

|1> _ 2' ~ |0,mk,t;0,mk,4.;0,m_k,t; l,m_k,4.), ( A - 2 ) 

1,+^,-k) = |0,mk,t;l)"ik,4.;0,m_k,t;0,m_k,4.), 
|1, — — , — lc) = |0, m kit ;0, m k ,4.; 1, m_ k ,t; 0, m_ k ,4.). (A.3) 
The eigenvalues are 

£i - E\ Q - h\ k . (A.4) 

Note that Ak is the coefficient of the kinetic energy term of the electron 
- see 10. 

3. With two electrons and total spin ±1 the total momentum of the electrons 
has to vanish because of Pauli's exclusion principle. The eigenvectors are 

|2, +1,0) = |l,TOk,t;0, m k ,4.; l,m_k,t;0,m_k,4.); 

2,-1,0) = |0,mk,t;l,^k,4.;0,m_k,t;l,m_k,j.). (A.5) 

The corresponding energy value is 

e 2 = E^ - 2hX k + hrj. (A.6) 

4. With two electrons and total spin zero the interaction with the harmonic 
oscillators contributes. An additional label is needed to distinguish be- 
tween the two levels of each of the Jaynes-Cummings pairs. The eigenvec- 
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tors are 

|2, 0, k, — > = sin^|l,m kit ;0,mk4;0,m_ k ,t; l,m_k,4.) 

+ cos^|0,m k ,t + 1; l,TOk,j.; l,m_k,t;0) m -k,4. - 1} 

(A.7) 

|2,0,k, +) = cos^|l,rnk,t;0,mk,4.;0,m_k,t;l!™-k,4,) 

-sin^|0,m k ,t - 1; l,«*k,;; 1, m_k,t; 0, m_k,4. + 1} 

(A.8) 

|2, 0, — k, — ) = sin^ k "|0,mk,t;l ! TOk,4.;l) m -k,t;0 ! m-k,4.) 

+ cos0 k ~|l,TO kit + l;0,rnk4;0,m_k,t; l,™-k,; ~ 1) 

(A.9) 

|2, 0, — k, +) = cos6» k |0,m kit ; l,m k ,j.; l,m_k,t;0,m_k,4.) 

-sin 6^ |l,m k)t - 1; 0, m k> ±; 0, m_ k)t ; l,m_ kji + 1). 

(A.10) 

The angles 9^ are given by 

ff£ = arctan T ^ ^ + ^ + 4 ^) » ( A ' U ) 

like in the Jaynes-Cummings model [2(j. The energies for these states are: 



4 = E* HO - 2hK ± - ^/K + w_ k )2 + 4^, (A.12) 

where the plus sign holds for the first and the third, and the minus sign 
for the second and the fourth class. 
5. With three electrons there are again four different classes of states: 

|3, +i k) = |l,m k)t ;l,m k) 4.;l,m_ kit ;0,m_ k ^) (A.13) 

|3, +i — k) = |l,m k . t ;0, m kii ; l,m_ k:t ; l,m_ k ,j.) (A. 14) 

|3, — — , — k) = |0,m k , t ;l,m ki ;;l,m_ k , t ;l,m_ k 4) (A.15) 

|3, — i k) = |l,m k:t ; l,m kii ;0, m_ k:t ; l,m_ kji ) (A.16) 
all with the same energy 



^3 = E* - 3/iA k . (A.17) 

6. With four electrons: 



|4, 0, 0) = |1, m k , t ; 1, m Kl ; 1, TO_ k , t ; 1, m_ k ,;}. (A. 18) 

The energy is 

e 4 = e)j - 4flA k + (A.19) 
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